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SECTION  I 


INTRODUCTION 


A  generalised  schematic  of  a  jet  engine  test  cell  is  shown 
in  Figure  1.  Generally,  the  width  of  a  plume  issuing  from  a  test 
cell  is  approximated  by  the  stack  dimension.  The  stacks  are 
usually  square;  however,  they  contain  acoustical  baffles  which 
considerably  reduce  the  open  area.  For  the  calculations  made 
in  this  report,  the  plume  width  was  assumed  to  be  the  square 
root  of  the  net  open  stack  area. 

Dimensions  can  vary  considerably,  depending  upon  the 
particular  cell  design,  but  the  principle  of  operation  is 
always  the  same.  An  engine  that  has  been  repaired,  or  otherwise 
maintained,  is  placed  in  the  cell  to  test  it  under  flight 
conditions  before  being  remounted  on  the  aircraft.  The  engine 
is  considered  a  mobile  emission  source  vhich  is  governed  by 
Federal  rather  than  state  regulations.  However,  the  test  cell 
is  immobile,  and  on  that  basis  a  U.S.  District  Court  upheld  the 
right  of  the  State  of  California  to  regulate  test  cell  emissions 
(Reference  1)  which  occasionally  violate  state  visibility 
requirements  of  Ringelmann  1  (20%  opacity) .  Since  the  U.S.  Air 
Force  has  a  large  number  of  test  cells  in  California,  this  court 
ruling  can  have  a  significant  impact  on  Air  Force  operations  and 
capital  expenditures. 

In  order  to  satisfy  state  regulations,  there  are  three 
possible  alternatives:  (1)  design  smokeless  engines  and  in¬ 
stall  them  on  all  existing  aircraft;  (2)  introduce  fuel  addi¬ 
tives  to  minimise  soot  formation;  (3)  use  particulate  control 
devices  to  treat  the  test  cell  exhaust. 

The  first  of  these  alternatives  is  already  being  pursued 
as  a  result  of  the  military  incentive  to  reduce  visibility  of 
in-flight  aircraft.  However,  replacement  is  very  costly  and 
time  consuming,  due  to  the  variety  and  number  of  existing 
aircraft  and  aircraft  engines. 

The  second  alternative  can  be  effective.  However,  fuel 
additives  are  organo-metallic  compounds  (e.g.,  Ferrocene), 
which  deposi.  metallic  oxides  on  engine  surfaces.  Considering 
the  cost  of  the  engine  and  its  maintenance,  and  the  cost  of  the 
aircraft,  anything  that  may  permanently  alter  engine  parts  is 
considered  highly  undesirable. 

The  third  alternative  does  not  affect  the  engine.  Because 
of  this,  it  is  the  only  alternative  which  state  regulatory 
authorities  can  impose.  Nevertheless,  particulate  control 
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Figure  1.  Generalised  Test  Cell  Schesutic 
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devices  cannot  be  designed  to  meet  an  opacity  requirement#  only 
a  specified  degree  of  particulate  removal. 

The  purpose  of  this  study  was  to  establish  the  connection 
between  test  cell  particulate  emissions  and  plume  visibility  as 
a  basis  for  specifying  control  devices  that  could  be  mounted  on 
the  test  cell  exhaust  stack.  Zn  addition#  theoretical  calcula¬ 
tions  were  made  to  see  under  what  conditions  electrostatic 
precipitators  or  venturi  scrubbers  might  satisfy  opacity  regu¬ 
lations. 
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SECTION  XX 


SMOKE  NUMBER 


Much  data  on  jat  angina  particulata  amissions  ara  in  tha 
form  of  S/VE  smoka  numbars  (SN) ,  which  maasura  tha  ralativa 
contrast  of  a  standard  filtar  papar  axposad  to  tha  axhaust 
amissions  for  a  standard  pariod  of  time.  A  faw  invastigators 
(Rafarancas  2-6)  hava  takan  simultanaous  maasuramants  of  par- 
ticla  loading  ("soot  dansity")  and  smoka  number.  Fewer  studies 
(Rafarancas  6-8)  have  datermined  plume  opacity  as  a  function  of 
smoke  number.  These  data*  which  ara  generally  of  poor  quality* 
arc*  plotted  in  Figures  2  and  3  for  various  anginas. 

Also  presented  in  Figures  2  and  3  arv~  empirical  cor¬ 
relations  based  on  Reference  9.  The  correlation  in  Figure  3 
includes  tha  results  of  Connor  and  Hodkinson's  work  (References 
10*  11)  relating  observed  visibility  to  plume  transmittance. 

As  can  be  seen,  there  is  about  as  much  error  in  predicting 
tha  Ringalmann  number  directly  from  the  smoke  number*  as  there 
is  in  predicting  mass  loading.  However*  in  order  to  determine 
opacity  from  loading*  the  particle  size  distribution  must  be 
known  (involving  an  additional  error)*  and  a  computer  program 
used  to  make  the  calculation.  In  either  case*  the  use  of  smoke 
numbers  is  a  very  unreliable  tool  in  predicting  test  cell  plume 
opacity. 


0  20  40  60  80 

SMS  Smoke  Number  (SN) 

Figure  2.  Relationship  Between  SAB  Smoke  Number  and 
root  Density  (Reference  6) 
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Correlating  Equations* 
ln(l  -  jfg)  -  - (.04681)®* ®}SN <40 

In  (1  -  j§5>  -  - (.03161) ®*®)SN*40 
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Ringelmann  Number  (RN) 


Figure  3.  SAE  Smoke  Number  vs.  Ringelmann  Reading  (Reference  6) 


Correlating  Equations: 


O  J79  (Reference  2} 

•  J79-GE-10  (Reference  8) 
Q  vT65  (clean) 

■  J57-P-8 

□  J57  (Reference  7) 

0 J52-P-6A 

A J52-P-6A  (Smokeless) 

A  7F30-P-8 

A  7F30-P-6  (Smokeless) 

#  756-A-10 


In(y^)  -  - (0.194L)  l-ln'l—j^)  )2*0®  SN  a  30 
(L  ■  plume  width  *  3.5  i) 

In^)  -  -(0.0617L)  l-lnd-^lf)),  SN  £  30 

RN  -  -.03757  +  4.5,  15  &  7  S  85 

RN  »  -.0800T  ♦  8.1,  95  2  7  >  85 

RN  ■  -.06677  +  5.0,  7  <  15 


SECTION  III 


VISIBILITY  EQUATIONS 


Jet  engine  exhaust  emissions  largely  consist  of  fine 
particles  of  unburned  carbon.  Because  they  are  black,  and 
therefore  absorb  much  of  the  incident  light  intensity,  carbon 
particles  will  exhibit  very  little  back-scattering  of  ambient 
light.  The  visibility  of  black  plumes  is  almost  entirely  a 
function  of  the  relative  contrast  between  the  background  sky¬ 
light  and  the  amount  of  such  lig» transmitted  through  the  plume 
‘  (Figure  4).  This  relative  contrast  is  independent  of  observer 
position  and  can  be  calculated  by  the  following  equation (Ref¬ 
erence  11) : 


T  •  ~ 


Bt 


exp 


where:  T  »  transmittance,  or  relative  plume  brightness 

*  brightness  of  light  transmitted  through  the 
plume 

B0  ■  background  sky  brightness 
W  ■  particle  loading 

R  *  plume  diameter 


P 


P  " 


I  - 


subscript  i 


particle  density 

total  number  of  particle  sizes 

particle  size 

extinction  coefficient,  or  ability  of  a  given 
particle  to  reduce  the  intensity  of  the  trans¬ 
mitted  light 

particle  size  in  the  distribution 


The  extinction  c* efficient  for  any  given  particle  is 
determined  by  the  following  equation  (References  11  and  12): 
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Black  Plumes 


For  transmitted  light*  the 
observer  sees  only  relative 
contrast  with  background  sky 
brightness o 


Figure  4. 


White  Plumes 

Scattered  light  originates 
from  both  in  front  and  be¬ 
hind  the  observer. 
Therefore*  the  relative 
position  of  the  sun  is 
significant. 


Visibility 
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Qext  m  JT  L  (2n  +  1) Real (an  +  bn) 

n 


where:  x  •  u  dp/X 

X  •  wavelength  of  light  in  which  the  plume  is  viewed 
(the  computer  program  allows  this  value  to  be 
input,  or  if  left  blank  assumes  an  average  value 
for  skylight  of  0.550  microns) 

an  &  b„  are  complex  Riccati-Bessel  function  of 
order  ''n": 


V <y)*n(x)  -  V (*> 

an  “  Vn'(y)tn(x)  -  mVj^y) Vn“*  (x) 


_  m4>n* (y)^n(x)  -  ^n(y)^n' (x) 
bn  “  ®4>n‘<yKn(x)  -  *B(y)?n‘l*) 


y  »  mx 

m  ■  complex  refractive  index.  This  is  a  function  of 
the  wavelength  of  light  (X)  at  which  it  is 
measured,  and  also  the  method  of  generating  the 
soot  particles.  The  computer  program  allows 
this  value  to  be  input,  or  if  left  blank  assumes 
a  value  for  amorphous  carbon  at  0.550  microns 
of:  1.96  -  0.66i  (References  13,  14).  (Note: 
Because  the  transmitted  light  is  altered  both 
in  phase  and  magnitude,  a  vector  is  needed  to 
express  the  effect  of  the  particles.  As  in 
electrical  engineering,  a  vector  can  be  expres¬ 
sed  as  a  complex  number  with  real  and  imaginary 
parts.  This  is  the  case  for  the  particle 
refractive  index.) 


Since  these  equations  involve  series  functions  of  complex 
numbers,  their  solution  is  not  simple.  instabilities  can 
easily  arise  (particularly  with  large  values  of  x) ,  which  cause 
the  extinction  coefficient  to  oscillate  wildly  and  even  produce 
negative  values.  In  order  to  generate  stable  functions,  the 
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following  method  is  used  (References  12,  15): 


d[ln<l>n(z)]  4»n'(z) 

Pn (z)  »  - - -  -  -1 - 

dz  ^(z) 

-  (I,  .  atlnCn(x)I  .  £„•(*) 

Qn<*>  *  - 3i  -£^y 


Therefore: 


(where  z  *  y  or  x) 

.  Vx>  rpn<y>  2  mPn<x)  , 

n  *n<x>  lPn(y)  -  mQn(x)J 


.  »n<x)  fmPn(y)  -  Pn(x), 

n  *  in(x)  lmPn(y)  -  Qn (x) 


Therefore: 


V  <*> 
V<x> 

PnU) 


-  Vi(«)  -  f  V*) 

■  f-n-l  (x*  ”  x 

.  Vi<*>  -  ?  v*> . 

ii»^Tzr 


Viu> 


n 

z 


Jv-l(z)  (v  -  1/2} 
Jv(z)  ~  t 

where  J  *  Bessel  function 
v  -  1/2  ■  n 

Using  Lentz's  continued  fraction  method^®). 

Jv  -1  1  Wj  | w; |  I *3 ,*2,*].  1 . . 

Jv  |  W2 1  I W3  .v»2 1 . . . . »  . . . . . 
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Convergence  is  reached  for  Jv-i/Jvwhen  |wp,....»W}i  in  the 
numerator  equals  |wp,....,W2l  i-n  the  denominator. 

Qn(x)  is  generated  by  the  following  recursion  formula: 


Qn<*> 


§  -  Qn-1 <x> 


(where  Q0(x) 


Since 


V-l<x>  .  Jv-l<x) 

"l^rnrr  jv(x>’ 


Finally, 


♦  n(*)  ■’J'n-l(x)  l (x)  J 

(where  »j>0(x)  ■  sin  x) 


5n(x)  ■  'i'ntx)  +  iBn(x) 


,  2n-l 


Bn<x)  »  (^^JSn-ltx)  -  Bn-2(x) 


(where  B0(x)  ■  cos  x,  and  Bj(x)  ■  +  sin  x) 


Following  Wiscombe  (Reference  12)  the  order  (n)  of  these 
functions  varies  from  1  to  N.  where: 


N  ■  x  +  4x1/3  +  2, 

N  -  x  +  4.05x3/3  +  2, 
N  «  x  +  4x3/3  +  1, 


These  equations  can  be  used  to  generate  the  requisite 
Riccati-Bessel  functions,  the  extinction  coefficient,  and  fi¬ 
nally  the  transmittance. 

Once  the  transmittance  (T)  is  calculated,  the  Ringelmann 
number  may  be  obtained  from  the. empirical  correlation  of  Connor 
and  Hodkinson  (References  10,  11),  Figure  5. 

A  computer  program  was  written  to  perform  these  calcu¬ 
lations.  The  FORTRAN  listing  and  sample  runs  are  given  in  the 
!  Appendices. 

t 

i 


y  >  4200 
8  <  x  <  4200 
x  <  8 
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Transmittance  {%) 


Figure 


5.  Black  Plume  Ringelmann  Number  Correlation 
with  Transmittance  (References  10,  11) 


SECTION  IV 


ELECTROSTATIC  PRECIPITATOR  (ESP)  EQUATIONS 


A  standard  mathematical  model  (Reference  16)  was  used  for 
predicting  the  ability  of  a  wire-and-plate  ESP  to  operate  as  a 
particle  control  device: 


wi 


8.85xl0"5  EcEpdpi  5 

- - -  (TT1: 


-AnW; K 

n.  =  1  -  exp (— V— ) 


where:  h  * 


m^  » 
Ap/Q  * 


wi  “ 

K  - 
K  - 


n  *  Eni«>i 

i 

overall  fractional  collection  efficiency 

fractional  collection  efficiency  for  particles 
of  size  dpi 

inlet  mass  fraction  of  particles  of  size  dpi 

specific  collection  area  of  the  ESP  (Ap  *  total 
collection  surface  in  m*?  and  Q  *  gas  flow  in 
m3/sec) 

theoretical  migration  velocity  of  particles  of 
size  dpi,  in  m/sec 

empirical  constant 

dielectric  constant  of  the  particles  (dimen¬ 
sionless) 


w  ■  gas  viscosity,  in  cp 


Ec  •  electric  field  strength  near  the  discharge 
electrodes,  in  kV/cm 


E-  •  electric  field  strength  near  the  collecting 
plates,  in  kV/cm 

dp  a  particles  size,  in  microns 


14 


This  model  is  based  on  a  field-charging  mechanism  and  is 
valid  for  particles  larger  thv.r.  0.5  microns. 

If  the  collection  efficiency  for  an  ESP  is  plotted  against 
the  particle  size,  the  resulting  curve  will  exhibit  a  minimum 
in  the  range  of  0.2  to  0.7  microns.  Above  that  range*  a  field¬ 
charging  mechanism  predominates  and  the  efficiency  declines 
with  particle  size.  Below  that  range,  a  diffusion-charging 
mechanism  predominates  and  the  efficiency  increases  with  de¬ 
creasing  particle  size.  Since  most  of  the  particles  emitted 
from  a  test  cell  are  smaller  than  0.2  urn,  a  large  K  value  of  600 
was  used  in  order  to  compensate  for  the  lack  of  a  diffusion¬ 
charging  mechanism  in  the  model  equations. 

In  order  to  reduce  the  amount  of  input  data  needed  to  run 
the  computer  program  and  avoid  a  ^rior  design  of  the  ESP,  the 
following  values  were  assumed. 


£  *  3,  for  carbon 

U  *  0.024  cp,  for  air  at  350°F  and  1  atm 

Ec  «  Ep  -  (4;5-  inches?  (2^54  cm/inch)  '  3*50  kV/cm 

(where  9  inches  is  ger.e:ally  used  as  the  plate  -  to-plate  spacing 
in  utility-type  ESP's,  with  a  secondary  voltage  of  40  kV) . 


The  computer  program  takes  the  uncontrolled  test  cell 
emission  data,  calculates  the  fractional  efficiencies,  and 
then  determines  tr.e  outlet  particle  size  distribution  and 
loading.  This  information  then  gees  to  the  visibility  portion 
of  the  program  where  the  outlet  Ringelmann  number  is  calcu¬ 
lated. 
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SECTION  V 


SCRUBBER  EQUATIONS 


A  standard  mathematical  model  (Reference  17)  was  used  to 
describe  the  particle  collection  efficiency  of  a  high  energy 
venturi  scrubber.  As  with  the  ESP  model,  this  also  required  an 
empirical  factor  (f)  to  make  the  model  agree  approximately  with 
actual  data: 


18  P2,  Ql  dft  2  i 

“  *i>  =  ''#(Pp)t^)(d~)  Cl  <<*if  +  °‘7> 


i  i  ,kif  4  °»7,  0.49 

"  1#4  ln  (  0.7  *  k*f  +  0.7 


)} 


This  equation  is  written  in  dimensionless  form,  and  therefore 
any  consistent  set  of  units  may  be  used: 


»  fractional  collection  efficiency  for  particles 
of  size  dpi 

Pj,  *  liquid  density 

Pp  ■  particle  density 

Ql  *  liquid  flow  rate 


Ci 


ki 


f 


gas  flow  rate 

Cunningham  correction  factor  for  gas  viscosity; 
for  particles  that  are  the  same  size  or  smaller 
than  the  mean  free  path  of  the  gas  molecules  (X) 

2  X  '°-44  dPi 

1  +  “■  {1.23  +  0.41  exp  ( - —  )]  (dimen- 

dPi  X 

sionless) 


Stoke3*  parameter 


cipPdPj2 

9ugdd 


_VgT  (dimensionless) 


0.5  (dimensionless  empirical  factor  based  on  the 
author’s  experience) 
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Vg<r  «  gas  velocity  in  the  venturi  throat 

A?T  *  pressure  drop  across  the  venturi  throat 

gc  ■  Newton's  Law  conversion  factor 

The  following  equations  require  specific  units: 

d<3  ■  mean  drop  size  in  the  venturi  throat  (Reference 
18} ,  in  microns 


°i  * 
Ojt  - 

* 

Ql  * 


1920  .FT.  , 

3-S9  ~ 

surface  tension  of  scrubbing  liquid ,  in  dynes/cm 
density  of  scrubbing  liquids  in  g/cc 
viscosity  of  scrubbing  liquid,  in  cp 


flow  rate  of  scrubbing  liquid,  in  gpm 


Qg  ■  gas  flow  rate,  in  cfm 
VgT  in  ft/sec 


X  ■  mean  free  path  of  gas  molecules,  in  microns 
3. 78  Wg 

V  ^gPg 

Ug  “  gas  viscosity,  in  cp 
Rg  ■  gas  pressure,  in  psia 
Pg  ■  gas  density,  in  lbm/ft^ 


Again,  in  order  to  minimize  the  input  data  requirements, 
the  following  operating  conditions  were  assumed: 

(1)  water  is  the  scrubbing  liquid  at  70°F 

(2)  gas  properties  are  those  of  air  at  350°F  and  1  atm 
pressure. 


For  venturi  scrubbers,  Qg,  /Qg  (the  1 iquid-to-gas  ratio)  is 
generally  5  to  30  gpm/1000  cfm,  and  £P<j<  is  10  to  70  inches  of 
water. 

The  computer  program  takes  the  uncontrolled  test  cell 
emission  data,  calculates  the  fractional  efficiencies,  and 
then  determines  the  outlet  particle  size  distribution  and 
loading.  This  information  then  goes  to  the  visibility  portion 
of  the  program  where  the  outlet  Ringelmann  number  is  calcu¬ 
lated. 

As  with  electrostatic  precipitators,  the  primary  collec¬ 
tion  mechanism  for  venturi  scrubbers  should  theoretically 
change  in  the  range  1.0  to  0.1  urn.  Above  1  urn,  the  particles' 
are  collected  by  an  inertial  mechanism,  while  below  0.1  urn  a 
diffusional  mechanism  should  prevail.  Again,  this  would  imply 
a  trough  in  the  fractional  efficiency  curve  for  particles  in  the 
1.0  to  0.1  pm  range.  However,  in  practice,  the  collection 
efficiency  of  venturi  scrubbers  continues  to  decline  below 
0.1  urn,  indicating  that  the  predominant  mechanism  remains 
inertial.  This  means  that  standard  venturi  scrubbers  are 
inherently  less  efficient  than  ESP's  in  collecting  particles 
smaller  than  0.5  ym.  One  method  of  overcoming  this  deficiency 
has  been  to  induce  condensation  in  the  gas  stream,  either  before 
the  scrubber  (by  quenching),  or  afterward  (by  utilizing  a  two- 
phase  ejector).  However,  these  methods  cannot  as  yet  be 
mathematically  modelled  with  any  confidence  and  have  not  been 
included  in  the  computer  program. 


SECTION  VI 


RESULTS  C.  DISCUSSION 


Grems  (Reference  19)  measured  the  particle  size  distri¬ 
bution,  loading,  and  transmittance  from  a  test  cell  at  Mc¬ 
Clellan  Air  Force  Base.  However,  the  particle  density  was 
unknown.  In  the  present  study,  this  density  was  used  as  an 
empirical  parameter  to  fit  the  computer  results  to  Greens'  data. 
Excellent  agreement  was  obtained  for  a  particle  density  of  0.92 
g/cc  (Table  1) .  For  comparison,  computer  results  are  also  shown 
for  *  particle  density  of  1.0  g/cc. 


Soot  particles  are  porous  spheres  of  carbon,  having  a  high 
void  fraction.  Solid  carbon  has  a  density  of  1.8  to  2.1  g/cc. 
Therefore,  a  particle  density  of  0.92  g/cc  implies  a  void 
fraction  of  about  0.53. 


Although  Grems'  data  indicated  a  bimodal  particle  size 
distribution,  a  straight-line  log-normal  fit  was  made  by  the 
computer  program  (see  Figure  6). 

Table  2  shows  computer  predictions  of  plume  visibility 
when  an  electrostatic  precipitator  or  venturi  scrubber  is  used. 
For  the  ESP,  a  specific  collection  surface  of  3281/m2  per  1000 
m3 /min  of  gas  (1000  ft2/1000  cfm)  corresponds  to  an  upper  limit 
in  commercial  applications.  Since  the  gas  flow  from  a  test  cell 
is  on  the  order  of  200,000  scfm,  200,000  ft2  of  collecting  plate 
would  be  required.  However,  even  such  a  large  ESP  has  only  a 
small  effect  on  the  Ringelmann  number  because  of  the  small 
particle  size.  The  size  range  which  has  the  greatest  effect  on 
visibility  matches  the  wavelength  of  visible  light— i.e.,  0.2  to 
0.7  microns.  This  is  precisely  the  range  in  which  an  ESP  or  any 
other  control  device,  is  least  efficient.  Therefore,  purchase 
of  an  electrostatic  precipitator  larger  than  the  test  cell 
would  only  have  a  marginal  effect  on  plume  visibility  during  the 
few  hours  per  week  the  test  cell  is  in  use. 

Similarly,  a  venturi  scrubber  designed  for  the  limit  of 
its  range  of  operability  would  also  have  a  marginal  effect  on 
test  cell  plume  visibility.  At  a  liquid-to-qas  ratio  of  30 
gpm/1000  cfm  (4.01  m^/min  water  per  1000  mVmin  gas)  the 
scrubber  would  produce  6000  gpm  of  waste  water  for  a  typical 
installation;  and  at  a  pressure  drop  of  70  inches  of  water  (131 
mm  mercury)  a  2400  hp  tail  fan  would  be  required.  Under  these 
conditions,  the  computer  predicts  an  overall  particle  collec¬ 
tion  efficiency  of  49%  with  a  visibility  improvement  of  0.5 
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TABLE  2.  COMPARISONS  OF  CONTROLLED  AND  UNCONTROLLED  EMISSIONS. 

UNCONTROLLED  EMISSIONS:  (Plume  Width  ■  8  Meter a) 


Particle  Loading 
(mg/m3) 

Particle 

Density 

(g/cc) 

Particle 

Size 

tun) 

Cum.Wt. 

%  Less 

Than 

Transmittance 

(%> 

Ringelmann 

Number 

6.34 

0.92 

21.9 

97.7 

66.2 

1.5-2. 5 

9.80 

97.7 

3.6S 

97.4 

* 

1.77 

95.8 

0.96 

94.8 

0.48 

94.2 

0.23 

88.2 

CONTROLLED  EMISSIONS t 

(1)  with  ESP:  SCA  «  3281  m2  /10C0  m3/nin 

Outlet  Particle  Collection  Transmittance  Ringelaann 
_ Loading  (mg/a>3)  Efficiency  (1) _ [%} _ Number 


2.93 


53.8 


82.8 


0.9-1. 9 


(2)  with  Venturi  Scrubber:  L/G  *  4.01  m3/min  water  per  1000  m3/min  Gas 

AP  *  131  mm  Bg 

Outlet  Particle  Collection  Transmittance  Ringelmann 
_ Loading  (mg/m3)  Efficiency  (%) _ _ Number _ 


3.26 


48.6 


81.0 


1. 0-2.0 
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Ringelmann  number.  Operating  data  with  a  scrubber  (not  a 
venturi)  at  the  Jacksonville  Naval  Air  Station  (Reference  20) 
indicate  an  average  particle  collection  efficiency  of  about 
75%.  However ,  there  was  considerable  uncertainty  in  the 
accuracy  of  the  data.  Stockham,  et  al.  (Reference  21)  also 
report  an  average  collection  efficiency  of  48%  with  water 
injection  into  the  augmentor  tube.  Again,  this  is  not  venturi 
scrubber  data,  but  it  does  indicate  the  validity  of  the  order 
of  magnitude  of  the  predicted  results.  It  should  be  noted  that 
opacity  measurement  with  a  scrubber  operating  is  virtually 
impossible,  since  the  scrubber  will  emit  a  large  and  obscuring 
steam  plume  of  its  own. 


/ 
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APPENDIX  A 


COMPUTER  INPUT  DATA  FORMAT 

The  following  pages  give  the  format  for  the  input  data 
needed  to  run  the  computer  program.  The  particle  loading  should 
be  determined  by  EPA  Method  5,  and  the  particle  size  distribu¬ 
tion  by  cascade  impactor  (Reference  19).  The  plume  width  can 
be  approximated  by  the  square  root  of  the  net  open  stack  area 
(the  actual  stack  cross-section  minus  the  area  occupied  by 
acoustical  baffles).  For  soot  particles,  the  refractive  index 
is  1.96-0.66i  at  a  wavelength  of  550  nm,  and  the  particle 
density  was  estimated  empirically  as  0.92  g/cc. 

Examples  are  given  for: 

(1)  Grams'  data^S) 

(2)  an  electrostatic  precipitator  with  a  specific  col¬ 
lection  area  of  1000  ft2/1000  cfm  (3281  m2  of 
collecting  plate  per  1000  m3/min  of  gas) 

(3)  a  venturi  scrubber  with  a  liquid-to-gas  ratio  of  30 
gpm/1000  cfm  (4.01  ra^/min  of  water  per  1000  mVmin 
of  gas)  at  a  pressure  drop  of  70  inches  of  water  (131 
am  Hg ) . 

CARD  »1 

This  card  contains  the  title  of  the  case  being  run, 
inserted  between  columns  9  and  10. 


CARD  #2 


This  card  contains  the  number  of  data  pairs  in  the 
particle  size  distribution.  The  minimum  number  is  2,  and  the 
maximum  is  100,  inserted  as  an  integer  between  columns  11  and 
15  (right  justified) . 

CARD  13  a,  b,*c,  etc. 

This  card(s)  contains  the  data  pairs  for  the  particle 
size  distribution.  Columns  1-10,  21-30,  41-50,  and  61-70 
contain  values  of  the  cumulative  weignt  percent  less  than 
;icle  size  dp;  while  columns  11-20,  31-40,  51-60,  and  71-80 
co  “ain  the  corresponding  values  of  dp.  Therefore,  a  maximum 
cf  four  data  pairs  can  fit  on  one  card.  If  there  are  more  data 
pairs  (as  per  CARD  #2) ,  these  are  put  on  subsequent  cards,  until 
the  tal  number  of  data  pairs  (cumulative  weight  percent  less 
than  dp,  and  dp)  is  equal  to  the  number  specified  in  CARD  #2. 
All  values  are  floating  point  numbers,  with  four  digits  (or 
blanks)  to  the  right  of  the  decimal  point  (right  justified) . 

CARD  #4 

This  card  contains  the  following  physical  parameters: 
columns  11-15 — the  effective  stack  diameter  in  meters 
(based  on  the  net  open  area)  expressed  as  a  floating 
point  number,  with  two  digits  (or  blanks)  to  the  right 
of  the  decimal  point  (right  justified) 
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columns  16-22 — the  particle  loading  in  mg/m3,  as  a  float¬ 
ing  point  number,  with  two  digits  (or  blanks)  to  the 
right  of  the  decimal  point  (right  justified) 
columns  23-28 — the  particle  density  ing/cm^,  as  a  float¬ 
ing  point  number,  with  two  digits  (or  blanks)  to  the 
right  of  the  decimal  point  (right  justified) 
columns  29-46 — the  particle  refractive  index.  Columns 
29-37  contain  the  real  part,  and  38-46  the  imaginary 
part;  both  as  floating  point  numbers  with  two  digits 
(or  blanks)  to  the  right  of  the  decimal  point  (right 
justified) .  If  these  columns  are  left  completely 
blank,  a  value  for  amorphous  carbon  of  1.96-0.66i  is 
assumed  by  the  program.  Note  that  the  refractive 
index  and  the  wavelength  that  follows  must  be  con¬ 
sistent 

columns  47-53 — the  wavelength  of  light  (in  microns)  at 
which  the  refractive  index  was  measured,  and  at  which 
the  plume  is  presumed  to  be  viewed.  This  must  be 
expressed  as  a  floating  point  number  with  three  digits 
(or  blanks)  to  the  right  of  the  decimal  point  (right 
justified)  .  If  columns  29-46  were  left  blank,  these 
columns  should  also  be  left  completely  blank,  in  which 
case  the  program  assumes  a  value  of  0.550  microns. 


CARD  #5 

This  card  contains  (in  column  2)  an  integer  number  which 
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indicates  whether  or  not  a  particulate  control  device  (elec¬ 
trostatic  precipitator  or  venturi  scrubber)  has  been  installed 
on  the  test  cell  exhaust: 

zero  (0)  means  no  control  device 

1  means  an  electrostatic  precipitator 

2  means  a  venturi  scrubber 
These  are  the  only  permissible  cases. 

CARD  »6 

This  card  depends  on  the  code  given  in  CARD  t5. 

(a)  If  there  is  no  control  device  (0  in  column  2  of  CARD 
#5) ,  CARD  f6  does  not  exist. 

(b)  If  an  electrostatic  precipitator  is  indicated  by 
CARD  *5,  CARD  #6  must  contain  the  specific  col¬ 
lection  area  (in  m2  of  plate  area  per  1000  a3 /min  of 
exhaust  gas)  in  columns  11**17  as  a  floating  point 
number  with  one  digit  (or  blank)  to  the  right  of  the 
decimal  point  (right  justified). 

(c)  If  a  venturi  scrubber  is  indicated  by  CARD  #5,  CARD 
#6  must  contain  the  liguid-to-gas  ratio  (in  m3/min 
water  per  1000  m3/min  exhaust  gas)  in  columns  11-15; 
and  the  scrubber  pressure  drop  (in  mm  mercury)  in 
columns  16-23;  both  expressed  as  floating  point 
numbers,  with  two  digits  (or  blanks)  to  the  right  of 
the  decimal  point  (right  justified) . 
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CARD  #7 


i 

l 

t 

;< ! 


This  card  contains  (in  column  2}  a  code  which  tells  the 
computer  if  more  cases  are  to  follow* 
sero  (0)  signifies  no  more  cases 
1  means  an  additional  case  follows 
For  each  additional  case,  CARDS  #1  to  7  must  be  repeated,  even 
if  some  of  the  data  remain  the  same. 


Figure  A-l.  Computer  Input  Data  Format 


Figure  1-A.  Continued.  Computer  Input  Data  Format 


Figure  A-l.  Continued.  Computer  Input  Data  Format 


APPENDIX  B 


OUTPUT  FORHAT 


Example  outputs  are  given  on  the  next  few  pages  for  the 
input  data  shown  in  the  previous  section  (i.e.  Grems'  data,  an 
ESP,  and  a  venturi  scrubber). 
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Figure  B— 1.  Continued .  Output  Format 


Figure  B-l.  continued.  CXitput  Fonaat 


APPENDIX  C 


FORTRAN  LISTING 

The  following  FORTRAN  program  performs  the  various  visi¬ 
bility  and  control  calculations.  It  was  originally  written  on 
a  UNIVAC  90/80-3,  and  later  converted  for  use  on  a  CDC  6600.  The 
version  shown  here  is  for  the  CDC  6600. 
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COMPUTER  PROGRAM  FLOW  CHART 


Figure  C-l. 
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READ 

SIZE 

DISTR. 

DATA 

CALL  GRAF 
(outlet  log¬ 
normal  size 

pis?*;. ) _ 


CALL  STACK 
(to  calc,  plume 
visibility) 


IF 
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Figure  C-2,  Computer  Source  Code  Listing 
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Figures  C-4,  Continued.  Computer  Soruce  Code  Listing 


Figure  C-2.  Continued.  Computer  Soruce  Code  Listing 
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Figure  C-2.  Continued.  Computer  Source  Code  Listing 


Figure  C-2.  Continued.  Computer  Source  Code  Listing 
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Figure  0-2.  Continued-  Computer  Source  Code  Listing 
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Figure  C-2.  Continued.  Computer  Source  Code  Listing 


APPENDIX  D 


Al 

A2 

A3 

AN 

AK 

AKF 

ALG 

B 

BIMAG 

B1 

B2 

B3 

BJ 

BJC 

BN 

CC 

CLCW 

CN 

CWDE 

CUDS 

CVJRD 


COMPUTER  PROGRAM  NOMENCLATURE 

real  part  of  an  denominator 
an  numerator /Al 
an  denominator/Al 

an 

Stokes'  parameter  (k)  in  Calvert's  scrubber  equa¬ 
tions  (dimensionless) 

AK (f )  +  0.7,  where  f  -  0.5 

liquid/gas  ratio 

intercept  of  straight-line  log-normal  equation  (in 
GRAF) 

i 

real  part  of  bn  denominator 
bn  numerator /B1 
bn  denominator /B1 
Jv-i/Jv  (with  real  argument) 

Jv-l/^v  (with  complex  argument) 

»>n 

Cunningham  correction  factor  in  scrubber  equations 
(dimensionless) 

cum.  wt.  %  in  1%  increments 

n/x 

outlet  cumulative  wt.  %  from  BSP 
outlet  cummulative  wt.  t  from  scrubber 
raw  cum.  wt.  %  input  data 
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DO 

■ 

scrubber  drop  size  (by  Nukiyama-Tanasawa 
tion) *  in  microns 

equa- 

DHUN 

is  defined  by  line  7  of  function  EINV 

DP 

m 

normalized  particle  size  from  subroutine  GRAF  (in 
microns) 

DPE 

m 

outlet  particle  size  from  ESP  (in  microns) 

DPRD 

m 

raw  particle  size  data  (in  microns) 

DPS 

m 

outlec  particle  size  from  scrubber  (in  microns) 

DPI 

m 

average  value  of  DP  in  interval  (microns) 

EFF 

m 

calculated  ESP  efficiency  (%) 

EFFS 

m 

calculated  scrubber  efficiency  (%) 

EINV(RCW) 

m 

inverse  normal  distribution  function  (RCW  • 

vari- 

able) 

is  defined  by  line  17  of  function  EINV 
is  defined  by  line  16  of  EINV 
-  1/3 

is  defined  by  line  15  of  function  EINV 
■  FLOAT (NVRD) 

is  defined  by  line  37  of  subroutine  ESP 

is  a  subroutine  that  uses  the  method  of  least 
squares  to  fit  the  raw  particle  size  distribution 
data  to  a  straight-line  equation 

I,IJ,IK,IL  are  all  counters 

ICASE  *  0,  when  there  are  no  further  cases*  and  ■  1  when 

another  set  of  data  cards  (i.e.  another  case) 
follows 

•  1  if  the  electrostatic  precipitator*  or  scrubber* 

has  such  a  high  efficiency  that  virtually  all  of 
the  original  particle  size  distribution  is  col¬ 
lected.  Otherwise*  it  is  0. 


ETA 

ETASQ 

E3 

EF2 

FND 

G 

GRAF 


INOGO 


IP 

■  p  (in  generating  values  of  w) 

IQtilP 

•  0,  if  there  is  no  particle  collection  device; 

if  there  is  an  electrostatic  precipitator;  - 
there  is  a  venturi  scrubber 

■  1. 
2,  if 

ITTLE 

•  alphanumeric  variable  (maximum  of  10  letters) 
responding  to  the  case  title. 

cor- 

JfJI 

are  counters 

K 

■  counter 

L 

■  (-1)P+1  (in  generating  values  of  w) 

N 

■  n  (order  of  Riccati-Bessel  functions) 

NE 

*  number  of  data  pairs  for  outlet  distribution 
ESP 

from 

NHAX 

-  max.  value  of  n 

NS 

*  number  of  data  pairs  for  outlet  distribution 
scrubber 

from 

NVRD 

*  number  of  raw  data  pairs  for  inlet  size  distribu¬ 
tion 

PATH 

*  plume  width  ■  effective  stack  diameter  (in  meters) 

PD 

■  pressure  drop  across  scrubber  (in  am.  mercury) 

PNX 

•  Pn(x)  *  ♦n' (x)/*n(x) 

PNY 

-  Pn(y)  -  V (y)/*n(y) 

PTLN 

■  natural  log  of  scrubber  penetration 

Q 

■  Qn<*)  “  *n'<x)/5n(x) 

QEXT 

■  extinction  coefficient  ■  Z(2n+1)  real  (an+bn) 

QSUM 

■  £(2n+l) real (an+bn) 

RB1X 

•  *nt«) 

RB2 

•  8n<*> 

RB3 

»  *n<*> 
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REFRAC 


RHOP 


RUM 


RIRL 


RNMAX 


RNMIN 


SLOPE 


STACK 

STEPB 


SOME 


*  complex  refractive  index  of  particle  (dimension¬ 
less) 

*  particle  density  (in  9/cc) 

*  imaginary  part  of  refractive  index 

*  real  part  of  refractive  index 
«  calculated  Ringelmann  number 

*  upper  estimate  of  Ringelmann  number  (■RNM+0.5) 

*  lower  estimate  of  RNM 

*  fractional  penetration 

0  specific  collection  area  of  ESP  (in  m*  per  1000 
m3 /min  of  gas) 

is  defined  by  lines  11  and  14  of  EINV 

*  slope  of  straight-line  log-normal  equation  (in 
GRAF) 

is  the  subroutine  that  calculates  plume  visibility 

-  1  (%SV 

i  dP  i 
**  «  penetration 

*  cumulative  wt  t  penetration  for  particles  of  sise 
DPE 


SUMS 


cumulative  wt  %  penetration  for  particles  of  size 
DPS 


SUMX 


SUMXY 


SUMX  2 


SUMY 


uxh 


TRAMS 


calculated  fractional  transmittance  ■ 

,.B[  f-3WD.  1  10r°  ,Qext.  , 

•xpnTpp>TOT  i-I1 
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TRANSW  *  %  transmittance 


VGT 

m 

gas  velocity  in  venturi  scrubber  throat 
ter s/sec) 

WD 

m 

|wp, . . . . ,*2l *  with  real  arguments 

NDC 

m 

|wp, . . . . ,W2l ,  with  complex  arguments 

WI 

M 

IMAG (WNC-WDC) 

WN 

m 

1 wp #••••» w^l#  with  real  arguments 

WNC 

• 

I  Wp, . . . .  ,W]_  I  #  with  complex  arguments 

WP 

at 

wp  -  (-1)P+1 

WPC 

m 

(-1)P+1 

WR 

9 

real  (WNC-WDC) 

X 

m 

EIMV(CWRD) 

XI 

m 

EINV (CLCW) 

X2 

m 

X2 

XLAM 

» 

wavelength  of  light  used  to  view  the  plume 
microns) 

XV 

m 

"dp/1  (dimensionless) 

y 

m 

logio  DPRD 

YV 

m 

m"dp/l  (dimensionless) 

ZI 

m 

particle  loading  (in  mg/m2) 
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